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1. Introduction 

It is Icnown tliat for two-dimensional Yang-Mills theory, the action TvF^ does not 
enjoy the uniqueness of its higher dimensional counterpart. In Migdal's lattice formulation 
, the local Boltzmann weight for a plaquette A of area a^ is a heat kernel for the operator 
TrE^, given by a sum over all irreducible representations of the gauge group: 

Za{U) = J2idimR)e-^"'^"'^^''^ Xr{U), (1.1) 

R 

where U is the holonomy of the gauge field around the plaquette, g is the gauge coupling 
and C{R) is the second Casimir of the representation R. Expanding about U = 1 one 
shows that this approximates the continuum Yang-Mills action arbitrarily well as a — *> 0, 
and because of the well-known self-reproducing property, even with finite area plaquettes 
the theory is exactly equivalent to a continuum theory. 

The self-reproducing property is true for an arbitrary function C{R). Such a lattice 
action ( p. . 1|) will correspond to a first-order continuum action in D = 2 of the form: 

S = j d^x [T,{E e^^F^,) - g^f{E)\ (1.2) 

Here F^i, = d^A^, — d^A^ + i[Af^, A^] is the Yang-Mills field strength and E is a, scalar field 
in the adjoint representation of the gauge group. Standard dimensional analysis applied to 
( p..2|) gives Ffj_i, dimension 2 and E dimension D — 2, so for D = 2 power counting allows an 
arbitrary function /. As has been pointed out 0, the self-reproducing property of ( |1.1| ) or 



the area-preserving diffeomorphism invariance of (|1 . 2|) allows exact solutions for arbitrary 
functions C and /, respectively. 

In this paper, we study the generalized two-dimensional SU{N) gauge theory i \L% ) 
minimally coupled to fundamental fermions in the large A^ limit. We demonstrate the 
consistency of this generalization and present ample evidence that it is both highly non- 
trivial and very natural. We thus add further examples to the vast number of known 
universality-classes of two-dimensional field theories and start addressing an amusing math- 
ematical question: How does this modification of the gauge interaction affect the natural 
observables in QCD2, namely the spectrum of mesons and their S'-matrix? But the main 
motivations for our work are two-fold: First, similar modifications of the bare QCD4 ac- 
tion might ultimately make this theory more tractable. Second, it is extremely important 
to keep looking for exactly solvable models whose analytical structure can be more fully 
understood or which more closely resemble four-dimensional QCD. 

Indeed, at present the analytically tractable model that is most analogous to QCD4 
seems to be QCD2 in the large N limit. |]1|] To support this point of view, we first point 
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out that no analytic description exist of the spectrum of D > 2 strongly coupled field 
theories. The situation is better in two dimensions. For example, the D = 2 sigma 
models with SU{N) target spaces have very interesting analogies with QCD4, notably 
logarithmic scaling, a mass gap and dimensional transmutation, however the spectrum 
and exact S-matrices show fewer analogies, and there is no direct analog of confinement. 
Another interesting possibility is to consider adding adjoint matter in QCD2, which gives 
a theory with roughly the number of degrees of freedom expected in a higher dimensional 
large A^ gauge theory [Q and may in the near future be analytically tractable. Finally, 
the finite N model is perhaps not much harder to work with numerically than the large 
N theory; adopting light-cone gauge and light-front quantization works equally magical 
simplifications at finite N, the essential difference is the presence of quark loops, which 
have been verified to be unimportant for certain interesting quantities (e.g. low lying 
masses). Despite this alternative, our own attitude at present is that we pay little for 
the additional simplification of large N. Furthermore we find the hopes for rewriting the 
theory as a string theory (which only makes real sense for large N) attractive. 

Pure gauge theory ( p.. 21 ) is exactly solvable for any N and has no local degrees of free- 
dom. When fermions are introduced, the non-quadratic nature of the function f{E) will 
introduce complicated interactions and a general model will no longer be solvable. Never- 
theless, since the only dimensionful coupling constant in the theory, g^, has positive dimen- 
sion, one expects that the generalized theory remains super-renormalizable. As an example, 
we will begin in section 2 by considering the simpler model of two-dimensional U{1) gauge 
theory with fermions but with more general gauge interactions as in ( |1 . 2| ) . This "general- 
ized Schwinger model" can be studied by the standard bosonization method. While the 
original (massless) Schwinger model effectively describes the meson (fermion-anti-fermion 
bound state) as a free massive scalar, the generalized model is an interacting scalar theory 
with a general Landau-Ginzburg potential. In two dimensions, Laudau-Ginzburg scalar 
theory is super-renormalizable, and leads to a rich family of critical universality classes, 
described by c < 1 conformal field theories. 

We continue our discussion in section 3 with a special class of generalized QCD2 in 



large N limit, with gauge field action (|1.2|) and f{E) = ^^ fnTrE'^. With a single Tr in 
( p..2|) , we are able to follow 't Hooft's original solution and derive the modified bound state 
equation as an explicit function of /. The infrared problem already present in the original 
model turns out to be even more subtle in our case and we have to carefully evaluate a 
class of highly singular loop integrals. After dealing with this technical problem we find 
that all infrared divergences cancel out of the bound state equation, thus giving strong 
evidence for the consistency of our models. The resulting modified equation turns out 



to be mathematically very natural. So far we are unable to present new exact solutions; 
however, we prove that the massless state present for the standard theory in the chiral 
limit (massless bare quarks) continues to be a solution for any function /, as is expected 
on general grounds. After discussing qualitatively how the linear force law is affected by 
the new terms in the action we go on to present the results of a preliminary numerical 
study of the simplest possible modification of the bound state equation, produced by taking 
/s 7^ 0. For small coupling the spectrum is a simple modification of that for the original 
theory and perturbatively (in 1/A^) the theory appears sensible and stable. The lowest 
massive meson state (present for massive bare quarks) is driven towards mass zero, and a 
critical coupling /3c exists with a new massless state. (This is in the free, leading order 
in 1/A^ theory - more generally this demonstrates the existence of critical points in these 
theories). For larger coupling the spectrum contains a tachyon. 

The potential relevance of our study to more physical dimensions D > 2 is discussed 
in section 4. If QCD2 is studied as a simple analog of QCD4, any of the generalized 
theories introduced above has a priori an equal right to be considered. More precisely, if 
an analytic technique exists which makes sense in arbitrary dimension 2 < D < 4 and 
whose continuum limit is unique in higher dimensions, then the dimensional continuation 
of this unique theory should give a preferred theory in D = 2 which may or may not be 
described by the TrF-^ action. 

Turning this around, if we find a member of the general class of two-dimensional 
theories with some special simplicity, perhaps this simplification has an analog in higher 
dimensions as well. We discuss the possibility of using the generalized actions in higher di- 
mensions, which generically produces a theory with the same continuum limit, but possibly 
more convenient cutoff scale dynamics. 

In particular we propose that to learn more about the connection with string theory 
it may be necessary to consider the general gauge theories, since a simple string theory 
may give a complicated spacetime theory. We will argue that the string picture that 
emerges from recent work on large A^ QCD2 [^^ can be simplified considerably if one 
allows generalized gauge interactions in the target spacetime. 

Section 5 contains concluding remarks. 

2. The generalized Schwinger model. 

The Schwinger model describes standard U{1) gauge theory minimally coupled to 
Dirac fermions in two space-time dimensions. While being exactly solvable (for the massless 
case), it demonstrates a rich spectrum of phenomena, such as the vacuum angle, quark 



confinement and the Higgs mechanism, and has been studied as a toy model for more 
physical theories in four dimensions. Detailed accounts of the model and its properties 
have been given in many different approaches H] and will not be reviewed here. We will 



instead consider its generalization in the sense of (|1.2| ), and point out some new features. 



The Lagrangian for the generalized Schwinger model is 

L = -Ee^'^F^, - g^f{E) + ^^{id^ - A^)^ - m^il;, (2.1) 

IT 

In this first-order formalism, the usual Maxwell theory corresponds to taking the scalar 
potential f{E) = 2E'^/n'^. A linear term in / would be the usual theta term. The field E 
is a pseudoscalar and thus if f{E) is not an even function, parity will be explicitly broken 
in the model. 

The field-theoretic aspects are most clearly exhibited by bosonization (originally dis- 
covered in this context). The standard rules are: 



1 

'^'0 = cm cos(2(/)) 



tl^iYd^tl;=—d^(Pd^<P, (2.2) 



where c is a constant related to the normal-ordering in defining the composite operator. 

In the second-order formalism of the usual Schwinger model, one chooses a convenient 
gauge (such as axial gauge) so that the gauge field equation contains no time derivative, 
and can be solved as a constraint equation. The same argument applies to the first-order 
formalism and a general potential f{E). The bosonized Lagrangian can be written as 

L = -e^^d^A, {E-4>)- g^f{E) + ^d^(t)d^4> - cm" cos(2(^) . (2.3) 

TT ZTT 

The gauge potential appears linearly and serves as a Lagrangian multiplier. Fixing a gauge 
Ai = and integrating out Aq simply gives a constraint that determines the scalar field 
E in terms of (j): 

E = <P+^- (2.4) 

where the constant 9 is the vacuum angle of the theory. After a convenient shift of (p by 
9/2, the final Lagrangian is: 

L = J-a^(/,a.(/) - g'^f{(j)) - cm^ cos(2(/) - 9) . (2.5) 

ZTT 

In the massless case (m = 0), the generalized Schwinger model simply describes an inter- 
acting scalar (the meson) with a general Laudau-Ginzburg potential /. 



From this example, we see that a generahzed two-dimensional gauge theory coupled 
with fermions remains super- renormalizable for a general potential /, as suggested by the 
positive dimension of the gauge coupling. The naive intuition that higher powers of F 
would be irrelevant is incorrect. The essential reason is that the gauge field has no local 
degrees of freedom. Power counting should be done with dimensions set by the fermion 
Lagrangian, giving A dimension 1 and E dimension 0. Unlike pure D = 2 Yang-Mills, we 
must give a renormalization prescription to completely define the model and determine the 
mapping between bare Lagrangians and physical models. A simple and standard choice 
would be to normal order the interaction using the free mass m bosonic propagator for 
contractions. 

Mesons exist as bound states of fermions, but they may have more complicated in- 
teractions determined by /. The qualitative physics of the generalized model is typically 
similar to the original model, but a number of interesting modifications are possible. 

First, there is the possibility of multiple vacua or even instability of the theory. The 
function f{E) gives the energy per unit length of a pure gauge field configuration of 
strength E. After bosonization, the function V{E) = f{E) + m? cos{2E — 6) has become 
a potential, and each minimum of V{E) is a possible vacuum of the theory. 

If we have multiple minima Ei, we have the possibility of stable states of non-zero 
charge q = rnr such that both left and right asymptotic fields are minima, if q = Ei — Ej 
for some i and j. The choice f{E) = cos2i? gives a particularly simple illustration. Each 
fermion now becomes a soliton with a gauge field 'dressing,' which modifies its mass while 
maintaining its non-interacting nature. 

Second, we can tune to critical points of the interacting scalar theory. The simplest 
critical point is attained by tuning the renormalized mass to zero in a theory with generic 
higher order interactions. Perturbatively this would even have been possible in the original 
massive Schwinger model by taking g^ negative; however such a theory would have energy 
unbounded below non-perturbatively. In the generalized model we can add a higher power 
of E with positive coefficient to fix this problem. The resulting theory has a phase transi- 
tion; for g^ < g1 the vacuum breaks the Z2 charge conjugation symmetry </> -^ —(p. Near 
g^ = g1 we see the critical behavior of the Ising model. 

The generalized Schwinger model should be a good qualitative quide to phenomena 
we can expect in generalized QCD2. A point to keep in mind however is that our large 
A^ treatment will only discuss the single meson sector. Clearly some of the field theoretic 
phenomena we described above will already have signals in this sector. For example, a 
critical point will be signaled by a meson mass going through zero as a function of the 
couplings. However, we will hardly be able to make a complete analysis just knowing the 



spectrum. The proper tool to study these phenomena would be the effective field theory 
for the mesons. 

Instabilities of the theory may or may not be visible in this sector. The reinterpretation 
of f{E) as a potential followed once the constant mode of 4> became dynamical, which is 
very much a field-theoretic effect. If we are constrained to the single meson sector, and 
have boundary conditions E = Eq for the the gauge field, we should only expect to see 
instabilities which can be detected by considering gauge fields with \E — Eo\ < q. This is a 
classical argument of course and the actual situation will be somewhat more complicated, 
as we will see. 



3. The generalized 't Hooft model. 

The Lagrangian is 

N N °° 

C = — TtE e^'^F^, - g^y^f^ TtE^ + i^ii^D^ - m)V^, (3.1) 

air 47r ^ — ' 

where E and the gauge potential A^ are N x N hermitian matrices. The field strength is 
F^i, = d^Aj, — d^A^ + i[A^, A^] and the covariant derivative is D^ = d^ + iA^. The A^'s 
have been introduced to make planar diagrams survive in the limit N ^ oo, with /„ and 
g fixed. The /^'s are dimensionless and g has dimensions of mass. We take one fiavor of 
fermionic quark in the fundamental. We have not taken the most general function f{E) of 
the introduction; powers of traces are also possible but will not be treated in this paper. 
If we take /2 = I/Stt, all other /n = and integrate out E we obtain the model of ^. 

As in m, taking light-cone gauge A_ = simplifies the theory tremendously. The 
preliminary observation that we eliminate the gauge field self-interaction is vital but this 
would also have been true in any axial or temporal gauge n ■ A = 0. What makes light- 
cone gauge simplest (and the gauge in which the theory has been solved most completely) 
is that we can then do light-front quantization of a theory with an instantaneous (in 
x^) interaction. In light- front quantization virtual pair creation is impossible, because of 
conservation of p_ and the positivity of p_ for every degree of freedom (assuming m > 0). 
Physical quark pair creation is subleading in 1/A^. Thus Hamiltonian evolution preserves 
the i/ji/j subsector of the Fock space and we can completely integrate out the gauge field - 
the quark self-energy is entirely reproduced by normal ordering. 

Although in some ways the light-front Hamiltonian description is more physical, the 
reduction to planar diagrams is clearer in Lagrangian perturbation theory, as used in |jl[], 
so our derivation of the bound-state integral equation will start there. 
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1. Planar Feynman rules in the light-cone gauge. 



In light-cone gauge, the Lagrangian reduces to * 



N 
An 



N 



TrEd.A+- — g^J2f-^^' 



n=2 



m 



(3.2) 



+V'Ja_^L + iA+(a+ + iA+)i^R - -T^ii^ti'R + i^Ri^L) 



and after solving for ?/'£,, 



- TrEd.A^ -Jza'Y. /-Tri?- + ^+(9+ + —- + iA+)^ 



Att 



n=2 



R- 



(3.3) 



We see that there are {EA) propagators, but no {EE) or {AA) propagator (we can 
regard Tri?^ as a vertex). Thus each propagator from an E'^ vertex must be connected to 
a quark line. This is what makes the model easily solvable - the gauge field self-interaction 
does not produce arbitrary fishnet diagrams but only a simple generalization of the rainbow 
and ladder diagrams of |]I[]. The Feynman rules are given in fig. 1, and a representative 
planar diagram is in fig. 2. 

We will discuss the qq bound state. (It is easy to see that the generalized interactions 
do not change the fact that exotic mesons are not bound in leading order in 1/N.) As 
for 't Hooft, the problem splits into two steps; evaluating the renormalized quark propa- 
gator, and then the renormalized Bethe-Salpeter kernel. Planar diagrams are generated 
recursively as in []I| , and the new elements in the quark self-energy are the "M" and higher 
order diagrams, while the 0{fn) correction to the kernel is a sum of graphs with I and 
n — I gauge legs attached to the quark and antiquark. All of these graphs can be expressed 
in terms of the "master integral" 
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absorb a factor of \/2i into ■0. 
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Fig. 2. A representative planar 



diagram. 



2 r°° 

TT J-oo 



^ -'5(A;i)- ^^5(^2)...^ —, 5(A;,-i) ^ 



P- - /Ci_ /Ci_ - /C2- A:n-2,- - /Sn-l,- /^n-l ,- - P'- 

(3.4) 
with n > 2 where S{p) is the renormahzed quark propagator 

S{P) = -. 2'^" v( ^^- ^3.5) 

Zp^P- — 171"^ — p-l [P-) + le 

The self-energy is 

CX) 

-iT{p) = -2/ 5^(-z)'^-in/./.(p,p) (3.6) 

n=2 

and the 2PI kernel is 



n-l 



K(g, q';p,p') = 4/ ^(-z)"7i/„ ^ /^(g, q')In-i{p' -q',p- q) (3.7) 



n=2 1=1 



where Ii{p,p') = l/{p — p'), and the incoming and outgoing quark (resp. antiquark) 
momenta are q and q' (resp. p — q and p' — q'). 

As in [0, it is important that the self-energy does not depend on p_|_ (as is clear from 
( ^.4| )). Physically this is because the interaction is instantaneous in a;"*". This allows us to 
do the /cj^-i- integrals independently without knowing F or m, producing 

/oo 
dki...dkn-i 
r 1 1 1 i^-*' 

-sgn ki- ;— sgn /c2 . . • sgn /Cn-r 



P-kl ki-k2 kn-2-kn-l kn-l-p' 



We have dropped the hght cone index { — } in (|3.8| ). This result could also be derived by 
normal ordering the terms of quadratic and higher order in the LC Hamiltonian. The 
sgn k- there comes from fermi statistics and the expansion of the field in modes with 
k- >0. 

The k- integrals are infrared divergent and must be regulated. The simplest prescrip- 
tion would be to follow Q and define the /c_ integrals as principal part integrals. However 
we will immediately run up against a subtlety first pointed out in this context by T. T. Wu 
T0[] : principal part integrals do not commute in general. A physical and safe procedure is 



to regulate the kernels -^-r by cutting them off at some small momentum A, as 't Hooft 
did originally [|l|. It is important to use the same X for all kernels in ( p78|) and safest to 
take A to zero only after performing all integrations.* A naive treatment of singularities 
by principal part integration in this problem leads to wrong results; i.e. they differ from 
the A — > limit of the finite A result. We found it convenient to invoke a slightly smoother 
version of 't Hooft's sharp cutoff; i.e. we regulate 

11/1 1 \ . s 

+ I—I^-T (3-9) 



p — p' 2\p — p' + iX p — p' — iX 
It is then straightforward to explicitly work out the first of the above regulated integrals 



I^ip,p')=l 


( 1 p-p' \ „p2 + A2 

K^p-p' {p-p')^+4Xy °^p'2 + A2 


— 


2A , p p\ 

1 i\i , ^x2(^^^*^^ \ +arctan 1 
(p-p') +4A2^ A A^ 


Taking now A to zero one finds 




Hp.p')-- 


2 V 

log TT^sgn p 5(p p) 
p — p p 



(3.10) 



(3.11) 

\j ■ 

It is not possible to express the regulated /^'s with n > 2 in terms of elementary functions. 
We must therefore find an alternative method to deduce the correct limit of these multiple 
singular integrals. It is convenient to introduce the generating function 

oo 

uip,p';z) = J2^~''UP:P') (3.12) 

n=0 



* It is by now clear that if properly derived from a physical regulator (which can be imple- 
mented in the original action), the principal part prescription, as used originally in [HJ^ and as 
will appear in our intermediate results below, is justified. One way to distinguish the situations 
in which it is valid is to note that in a Hilbert transform acting on continous functions of an 
appropriate type, orders of integration can be exchanged. In particular, the p.p. integral in our 
bound-state equation will act on smooth wave functions. A completely independent check of the 



spectrum found in ||I|] has also been done in |11], using the axial gauge treatment of [12|. 



and define 

(3.13) 



Io{p,p') -- 


=sgn p 5{p - 


1 
-P 


h{p,p) -- 


V 
p-p' 





where V denotes the Cauchy principal value. Formally, the generating function u{p,p';z) 
satisfies the integral equation 

z u(p,p';z) = z sgn p S(p — p') + / dk sgn /c u(k,p';z) (3-14) 

y_oo p-k 

We can give it meaning by noting that sgn p J dk^^ should be a well defined integral 
operator on continuous functions; in particular the functions sgn k \k\'^'^ are eigenfunctions 
of this operator: 

/•°° -p 

sgn p / dk -sgn k \k\'^'^ = n cot ttz/ sgn p Ipl^'' (3.15) 

i_oo P-k 

In view of eqs. ([3^ , ( |3.12| ) these functions are also eigenfunctions of the integral operator 
sgn p f dk u{j>, k; z): 

sgn p dk u{p, k; z) sgn k \k\^'^ = sgn p \p\^^ (3.16) 

J Z — 71 cot ITU 

We therefore have to find a kernel which reproduces (|3.16|) . This is done by noticing that 
sgn p j dk-^ I fcP" ^g'^ k \k\'^^ = TicotTi{y — a) sgn p \p\'^^ . Using the addition formula 
for cot7r(z^ — a) one finds the (unique) expression 

«(p,y;^) = ^^[^7 l4l'"^'^+^sgnp5(p-/)] (3.17) 

z^ + n"^ ^p — p' p' -' 

where z = ti cot Txa{z)^ and thus 2a{z) = — arctan -. One can check by direct computation, 
using the formula (derived from ( |3.9| )) 



^ ^ ^-{^ + T:^)-7v'd{p-p')l{d{p-k) + d{k-p')) (3.18) 



p — k k — p' p — pp — k k — p 2 

that u{p,p'; z) is indeed the solution of the integral equation (|3.14|) . Expanding in z~^ to 
C(2-4) one finds in addition to (P^),(pi]) 

h{p,P)=^ 7 log I — l-TT 



p-p' p' p-p' 

h{p,P) =- 7 log I — I - -TV 7 log I — I +7r sgnp d{p-p) 

6 p — p p 3 p-p p 
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It is important to notice that these loop integrals are generalized functions (distributions). 
One easily verifies to low order using these explicit formulae together with ( p.l8|) that we 
can recursively calculate the higher I^s from lower orders and that the result does not 
depend on the order of integration. In order to test the independence of our results on the 
details of the infrared cutoff procedure we numerically checked our formulae to this order 
using 't Hooft's sharp cutoff. We found convergence to the above result ( |3.19| ) as the sharp 
cutoff tends to zero. 

In order to calculate the contributions to the self-energy —ir{p) we take the limit p -^ 
p' of the distributions In{p,p')- AH terms with powers of log |^| higher than one tend to 
zero while the terms linear in the logarithm give a finite, non-zero 'mass renormalization'.* 



The terms -^-r tend to zero likewise due to our prescription 



0. We 



p—p' - - - -— - - jT jr - \| ^ y p—p' 'P~P 

conclude from ( |3.17| ) that the odd order mass renormalizations are zero. The even order 
self-energies are both finitely and infinitely renormalized; the latter due to the presence of 
terms ~ 5(0). Such infinite self-energies are already present in the original 't Hooft model 
and we have to study the full bound-state equation before getting worried. Following again 
Ij^ it is given by 

[2p+-m'{- + -^)](t){q)=[T{q)+T{p-q)]ct>{q)+ f dk K{q,k-p,p) (t>{k) (3.20) 
q p-q" Jo 

All momenta without index are to be understood as minus-components. We can express 
this equation through the generating function u{p,p'; z); this leads for the right-hand side 
to the result 



2 I dz ^,, 



^« * 2^^ '^' 



- {u{q, q; iz) + u{p~q,p- q; iz)) (t){q) + 



rP 
-f 2 / dk u{q, k; iz) u{p — k,p — q; iz) (f){k) 
Jo 



(3.21) 



where f{z) is as in (|3.1|) . Using the identity 

^q — k' {q — k)' 



+ 7r2 (5(0) 5{q-k) 



(3.22) 



* Strictly speaking, this is not a mass renormalization, as is clear from considering a zero 
bare mass quark, which would be protected from perturbative mass renormahzation by chiral 
symmetry. The choice of LC gauge breaks Lorentz invariance in gauge-variant quantities and 
allows producing a self-energy which enters the same way as a mass in the subsequent analysis. 
We thank N. Seiberg for a discussion on this point. 
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which is again derived from ( p.9| ) we can rewrite after some algebra the expression ( |3.21|) : 

f (z) -izy^ ,^2a{iz)[- ^ ) (p{q) 

i 7,^ — TT^]^ ^ n n — n' 

(3.23) 



2<7^ 



27X1 



[z^ 



- n^ 
2z^ 



(^2-7r2) 



2^2 



dk 



V 



Mp-k) .2oc{i 



{q — k)"^ {p — Q)k 



<P{k) 



AU the infrared divergences, i.e. aU terms with 5(0) have canceled! This is the generalization 
of the correponding phenomenon in the ordinary 't Hooft model. It serves as a highly 
nontrivial check on the internal consistency of our models. Note also that a further mass 
renormalization coming from the kernel is seen in ( p.23|) ; it remains however true that only 
the even couplings in the action lead to mass renormalization. After going to dimensionless 
variables in the usual way, 7 = tt^, ^u^ = ^2p_|_p_ and x = -, y = - one obtains the 
final form of the bound-state equation: 

^■)] 't'i^) = 



r 2 A 1 

'- ^X 1 



z^ + n^ 



iz- 



{z^ 



TV 



5p 2"(-) C- + ^) 



[X 



(3.24) 



2z^ 



71 



2^2 



(x-y)'^ 



x{l-y) 
(1 -x)y 



2a{iz) 



<t^{y) 



Eq. (|3.24| ) constitutes the principal technical result of the present work. The generalized 
kernel results in logarithmic corrections to the 't Hooft equation; for example, turning on 
the couplings /s and /4 in addition to /2 = I/Stt in the 't Hooft model one obtains 



ix^^{x) = (7 - 1 



160 



-TX 



dy- 



\x 



3 «(^r^)*w 



(3.25) 



327X U {log^^-^f] <Piy) 



x)y' 



Note the factor of i in front of the term linear in the logarithm, which came from the k^ 
integral producing (|3.4|) . It serves to ensure the hermiticity of our Hamiltonian. Indeed, 
as is easily seen from ( p.24|) , all terms with odd powers of the logarithm (corresponding to 
odd potentials) come with such a factor of i, while even powers are real. 

To complete the discussion we must find boundary conditions for (p{x) compatible 
with self-adjointness of the Hamiltonian implied by ( p.24| ) . As in the 't Hooft model, the 
appropriate boundary condition is (f){x) 



x^ at X = and 



[X 



1 — x)^ at X = 1, where 
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the exponent (3 is determined by requiring the leading term x^ ^ in the integral equation 
( p.24|) to cancel. It must satisfy 



/dz , z^ + 71^ 2z^ 

27n I (^^ — TT^)^ [Z"^ — TT^j^ 

(3.26 
Let us discuss some of the properties of the bound-state equation. A very useful 
identity is 

I'dy—^ {^'' = -2nu cot27r^ {- + j^) C-^)"" (3-27) 

Jq [x — yy y ^x 1 — x X 

It may be used to show with a little bit of algebra that for massless quarks (m^ = 0, i.e. 
7 = 0) the wavefunction (f){x) = 1 solves the bound-state equation with mass eigenvalue 
11^ = for an arbitrary potential f{z). We conclude that the massless ground-state is 
invariant under the generalized interactions! While being true in general, we can make 
this property manifest for a purely even potential by rewriting (|3.24|) for that subclass of 
models as 

dz^,,. 2z' f\ V rx(l-y)^'"('^) ^^-^^^ 



2^ <P 7T-n^) .... _.^. dy 



27Ti {z"^ — Tx'^Y Jo (^ ~ y) 



(1 -x)y 



((/>(y) - (/>(X)) 



The existence of such a massless state in the 't Hooft model is a consequence of 
chiral U{1) symmetry as shown in |jl^. Although in (|3.1| ) both SU{N) and U{1) vector 
symmetry are gauged, the large A^ limit takes the U{1) gauge coupling to zero in a way that 
suppresses the chiral U{1) anomaly. The simplest argument for this is that the anomaly 
receives contributions from Feynman diagrams with one fermion loop, which by standard 
large N counting is 0{1/N). This remains true in the generalized models. Conservation of 
the vector U{1) current in two dimensions allows writing J^ = e^'^d^cp, and conservation 
of the chiral U{1) current J^ = e^^J'^ then implies the existence of a free massless boson 

0- 



A more detailed discussion using non-abelian bosonization is given in [0. The 
fermions are represented by a /c = 1 SU{N) WZW model with a massless boson repre- 
senting the remaining U{1). The complete theory consisted of a minimally gauged SU{N) 
WZW model and a massless boson (p coupled with the U{1) gauge field TrF through (pTicF. 
For the 't Hooft model there is no coupling between SU{N) and U{1) sectors. In the gen- 
eralized models, terms TrE'^ with n > 3 produce explicit couplings between the SU{N) 
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and U{1) gauge fields at subleading order in 1/A^. This surely implies that at finite A^ 
long-distance physics eliminates the massless boson, and the bosonized action would be a 
good starting point for analyzing this. 

We have found that adding TtE"^ to the action adds new terms to the integral kernel 
in the bound-state equation up to -r-^-rw log"~ {p/p')- This is essentially the Fourier 
transform of the interquark potential. Only the even terms TtE'^'^ change the strength of 
the linear confining potential, while odd terms Tri?^""*"^ produce short range forces. The 
long-distance behavior of the potential can be understood classically. Solving Gauss' law 
d-E = jjJ- for a source at x" = produces E = qe{x~) (with e{x) = 9{x) — l/2). Solving 
the equation of motion d-Aj^ = f'{E) (with boundary conditions such that the background 
potential is zero) then reproduces qualitatively the behavior we found. However, the precise 
strength of the induced linear potential is not given correctly by such a simple argument 
for the new interactions n > 4. One reason for this is that it does not reproduce the planar 
nature of the large A^ perturbation theory. 

The new terms in the bound-state kernel are less singular at k = k' and correspond 
to short range forces. Clearly they can only be understood in the quantum theory. The 
terms .^^ ^ log"~ (x(l — y)/y{l — x)) with n > 2 in the bound-state kernel are analytic 
at a; — y = (zero momentum transfer) for x and y away from the endpoints and 1. If 
we neglected the endpoint singularities, their Fourier transforms to position space would 
therefore be contact terms. In the limit of heavy quark masses this would be justified; the 
wavefunction would be highly peaked (for equal quark masses) oXx = 1/2. The naive four- 
quark operators corresponding to these contact terms would be non-renormalizable even in 
two dimensions, but in the sense of an effective theory they would give a good description. 
Away from the heavy quark mass limit, the endpoint singularities become important, and 
it is not clear to us whether a local effective field theory description exists. Related to this, 
since the kernel is not only a function of momentum transfer, the interactions also have 
the slightly strange feature that the positions a;~ of the quarks can change instantaneously 
in x'^ during the interaction. 

The new interactions will shift the masses of all mesons except the massless one 
(present for vaq = 0). Using perturbation theory to estimate this, we find that for an 
even perturbation Tri?^"^ there is a correction at first order, with the same sign as the 
coupling, while for an odd perturbation TrE'^^'^^ the first correction is at second order 
and always negative. A WKB argument shows that for the highly excited states m^ ~ n 
the log coupling shifts rn^ by 0(?i^~^). We therefore expect that for sufficiently strong 
coupling (negative for even perturbations) we can drive a meson massless or even tachyonic. 
Classically one certainly expects instability for certain potentials, as we mentioned earlier. 
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We have preliminary numerical results for the case fs ^ 0, at the special mass 7 = 1. 
Our numerical method is based on the one developed in Hanson et. al. |jT^ and is described 
in the appendix. For the 't Hooft model with coupling /2 = I/Stt^ (chosen to eliminate 
some tt's from the following) and 7 = 1, the spectrum is (combining known analytic and 
numerical information) m^ ~ 71 — 0.25 — 0.01/n^ + . . . . 

Our results for /a = e/247r^ 7^ are consistent (at the 1% level for e < 0.3) with the 
slope remaining unchanged, and the intercept 0.75 decreasing to 0.75 — 1.04e^ for small e, 
consistent with second-order perturbation theory. There is a critical coupling Cc ~ 0.45 
where the lightest meson goes through zero mass. We interpret this as the instability we 
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discussed in section 2. For |e| < Cc we see no sign in this calculation of instability, though 
one might expect it to show up in a multi-meson sector. The main point is simply that 
the model in this sector is well-defined and different from the original 't Hooft model (for 
any quark mass). 

We expect models with f2n > 0, /2T1+1 = to be completely well-defined and stable. 
More generally, we expect that f{E) bounded below suffices, if we derive our bound state 
equation by expanding about the true minimum. 

For very large /s, it is not clear whether the model will confine. Standard techniques 
for singular integral equations |T^ should suffice to solve the pure /s model (since the kernel 



has a single pole) and one expects the integral operator in this case to have continuous 
spectrum. Very naively, a small /2 coupling looks like a singular perturbation which could 
restore the discrete spectrum and possibly confinement. However, it could be that the 



correct boundary condition here is given by a different branch of (|3.26|) than one takes for 
small /a, which might still allow continous spectrum (as one expects from the large mass 
limit.) 

As one approaches the critical value ec, the end-point behavior of the wave functions 
appears to approach x^ with purely imaginary (3. This leads to very poor convergence 
in the numerical treatment of the appendix and to check our results we repeated the 
calculation using a simple Taylor series basis (around x = 1/2) for the wavefunctions 
and Richardson extrapolation in the basis size. This produced consistent results for small 
e, but a much larger critical coupling e^ ~ 0.75. It seems likely that €c is determined 
by the condition Re /3 = and solving for this condition in (|3.26| ) numerically produces 
ec = 0.676215 ± 0.000005. 

As one decreases the quark mass, there will be an interesting and subtle crossover 
between the behavior we just saw, and the case of a massless meson, which is unaffected 
by these perturbations. It would be interesting to know if the resulting critical points are 
again roughly as expected from classical considerations. This requires more care in the 
niimerics and we reluctantly leave the question for future work. 



4. Higher dimensions 

It is conceivable that these generalized actions could be of direct value in higher 
dimensions. First, we hasten to reassure the reader that in D > 2, expanding around the 
standard (free) UV fixed point, the new terms Tri?" with n > 2 are non-renormalizable.* 



The only exception is the term ef,„xTrE''[E'' ,E^] in D = 3. 
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Rather than trying to define them in the continuum, we take as our definition of generahzed 
QCDx) a lattice gauge theory where for each plaquette Boltzmann weight we use ( |1 . Ij ) . 

The study of these generahzed theories in higher dimensions seems to go against two 
very widely accepted principles. The first is universality: the choice of bare action within 
wide limits is irrelevant to the study of the continuum theory, so why should we not be 
satisfied with the actions we know? We would reply with a two dimensional analogy: 
there are many models whose long distance limit realizes the physics of the three-state 
Potts model, for example, but if one chooses to study Baxter's definition, one can find 
exact solutions and make much quicker progress. Although one is not optimistic about 
exact solvability in realistic higher dimensional models, valuable technical simplifications 
might be possible. 

The second is that theories with higher derivative Lagrangians are non-unitary. We 
have already shown that this is not true in our generalized QCD2 but this might be thought 
to be a special case, unitary because the gauge field has no local degrees of freedom. 
However this is not the general explanation. One element of the continuum explanation 
is that one can write the corresponding continuum Lagrangians using higher powers of 
the field strength F^, but not higher derivatives of F, therefore preserving the number of 
canonical degrees of freedom. 

We now argue that generalized QCDd on a hypercubic lattice can satisfy reflection 
positivity and then is the Euclidean continuation of a unitary quantum field theory. As is 
well known it suffices to have a positive, self-adjoint transfer matrix. Choosing a lattice 
axis as time, the transfer matrix could be taken as the composition of an "electric" and 
"magnetic" operator, 

T{At,a)=TE{At,a)TB{a) (4.1) 

where a is the spatial lattice spacing and At the time step. (We will suppress the a 
dependence in the formulas below.) We will take this to act on gauge-invariant wave 
functionals of the spatial link variables. 

Clearly a more symmetric definition would be 

T{At) = TE{Atf''^TBTE{AtY/'^ (4.2) 

and it would suffice to show that T£;(At)^/^ and Tb are positive and self-adjoint. Since Tb 
is simply a multiplication operator by the product of the magnetic Boltzmann weights, we 
need the individual Boltzmann weights to be positive. As for TE'(At)^/^, in fact we can 
write 

TEiAtf" = TE{At/2) = exp -^ (4.3) 
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for a self-adjoint Hamiltonian He- This simply follows from the original definition of the 
generalized heat kernel action for the Hamiltonian H: we can show that 

^TE(t) = -J2H^TE{t) (4.4) 

i 

where Hi is the Hamiltonian for the site i. In detail: we have 

iTEm)iU)= f]JdUu]JdUl^]JZ,tiU,jU-^U;,Uu)ijiU'). (4.5) 

i ij ij 

Now the point is that 

^Zt{U V) = -H{Eu)Zt{U V) (4.6) 

at 

where the Hamiltonian is constructed from generators of left rotations of U , for all V . We 

can therefore evaluate the r.h.s. of ( ^.4| ) under the integral. 

We have shown that a generalized QCDo has reflection positivity (at finite cutoff a), if 
the solution of the generalized heat equation is positive for all U at time o?. The proviso is 
necessary for positivity of the magnetic plaquette Boltzmann weights and is not required 
in D = 2. That reflection positivity implies unitarity after continuation to Minkowski 
space-time is clearest in Aq = gauge, i.e. with the links Uu = 1. This choice would 
simplify ( |4.5| ) but have no effect on the argument. 

The condition of positivity of the solution of the generalized heat equation is quite non- 
trivial and in general is false. Typically, positivity will require conditions like tk << ^2" 
for k > 2. We have verified for the U{1) case that there are finite values of the couplings 
compatible with positivity and have no reason to doubt it for SU{N). Even if this positivity 
is lost, it is not obvious that the Minkowski theory is nonsensical, because we always have 
a self-adjoint transfer matrix. Correlation functions of operators whose time separations 
are even multiples of the lattice spacing will still satisfy reflection positivity. [^ Whether 
this observation is relevant probably depends on the details of a given case. 

A context in which a modified action has been proposed in the past is the string 
interpretation of the large A^ limit of the strong coupling expansion as derived in [^ . The 
motivation can be illustrated in D = 2 and although we allude to one could make 
the same point starting from the conceptually similar but simpler approach of Gross and 
Taylor 0. Without going into a detailed description of these results, the basic idea is to 
expand the plaquette Boltzmann weight as 

Z{U) = ^ JJ(A*Trt/*)'"^(l + 0(l/iV))exp-^2 j ^ |„.| + 0(i/iv) J (4.7) 

{rii} i \ i / 

18 



(by rearranging the character expansion) and interpret a term ]^(TrC7*)"^" as the contribu- 
tion of a string configuration which locaUy has rii world-sheets with boundary wrapping i 
times around the boundary of the plaquette, in other words a n-fold cover of the plaquette 
with n = ^^ \ni\. The leading term in the exponential then is interpreted as a dependence 
on the world-sheet area A as exp—g^A, while the 0{1/N) terms (which can contribute 
to leading orders of the free energy) are interpreted as due to insertions of additional fea- 
tures at branch points or other singularities of the covering. This is where the modified 
actions we are considering become relevant - changes to the action translate directly into 
changes to the 0{1/N) corrections of the action, and hence the additional world-sheet 
features. The idea is that from the point of view of this expansion, the simplest theory is 
not derived from the TrF^ action but instead from an action which eliminates the 0{1/N) 
corrections in the exponential. In the language of this corresponds to eliminating all 
world-sheet features except for the O and 0~^ points. This theory has been discussed (to 



a limited extent) in ^ and [|T8|. 



The most useful descriptions of the action are either as a particular case of ( |1 . 1| ) with 
f{R) = ^j \ni\, or as an expression similar to that of [jl^ in which the eigenvalues of U 
behave as free relativistic fermions and terms like exp— 0^/t are replaced by l/(^^ + 1^). 
In the form ( |1.2|) it apparently cannot be expanded around -E = 0, the usual starting point 
for perturbation theory. Nevertheless in the large N limit a number of calculations can be 
done with it and the results are non-singular. In the calculations which can be reduced to 
free fermions, the explanation is clearly that the fermi surface is at a 'momentum' E r^ N 
and we never see the E = Q singularity. However this is not a very general argument and 
one would like to see an analogous phenomenon in calculations using other methods. One 
of the original motivations for the present work was to examine this question; however 
the results presented here do not suffice, as this action requires products of the invariants 
Tr£^" for its expression. So far we see no reason why such models should be fundamentally 
harder to solve, but we leave such questions for future work. 

5. Conclusions. 

In this work we have formulated a generalization of two-dimensional gauge theory 
coupled to matter, argued that the generalized models are renormalizable and unitary, and 
using the techniques of the original work of 't Hooft, solved a large subclass of such models 
in the large A^ limit, in the sense that we derived a integral equation which determines 
the spectrum of mesons. The simple classical picture of the linear potential produced by 
a quark is qualitatively valid for the generalized interactions Tri?^"^, but all generalized 
interactions produce additional short-range forces. 
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The models are qualitatively similar to the original model, with a typical spectrum 
m? ~ n for sufficiently massive states, but can show new behaviors such as phase transitions 
or deconfinement. This was illustrated in a preliminary study of the model TrEF + TrE'^ + 
eTrE^ + i{j{p + m)'i/j. The results also have bearing on some old questions, such as the 
validity of the principal part prescription as an infrared regulator (in general, it is not). 

We foresee two types of applications for these models. The first is the possibility 
of choosing the action to produce a closer analogy to some feature of higher-dimensional 
gauge theory, such as logarithmic violation of scaling. The second, potentially quite im- 
portant application, would be if a model in this large class turned out to have qualitatively 
similar physics to the original model but was in some sense exactly solvable. Perhaps most 
interesting would be a model which was not integrable in the sense of having a factorized 
S- matrix yet which allowed analytic calculation of the S- matrix to any order. Presumably 
this could be done if closed expressions for the meson wavefunctions could be found. 

Light-front quantization is of considerable interest as a non-perturbative technique 
in higher dimensions, [^ and these models provide new toy examples involving massless 
states, symmetry breaking, and vacuum instability. 

Given the interconnectedness of physics and more specifically two-dimensional field 
theory, and the important role played by gauge fields, it seems safe to predict that unfore- 
seen applications will also be found. 

It is also conceivable that these generalizations have some value in higher dimensions. 
This is not the first time such a generalization has been proposed, but we believe the present 
work significantly clarifies the questions which would need to be answered to justify the 
use of such an action. 
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Appendix A. Numerical methods 



Our numerical method is essentially the one developed in |T^. We expand the wave 
function in a basis of Chebyshev polynomials: 

"-' (A.l) 

n 
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where Un{cos9) = sin(n + 1)9/ sin 9 and we have redefined x to make the hmits of inte- 
gration —1 and 1. The weight factor is not the one which makes the basis orthonormal 
but rather was chosen so that we can apply a standard formula from the theory of integral 
equations pO| : 

[ -^{1 - y^f/^U^-^{y) = -TrT^ix). (A.2) 

J y -X 

This weight factor also gives the prescribed endpoint behavior (3 = 1/2 for the special 
case 7 = 1. The basis is also quite suitable for (3 > 1/2, however it does a bad job at 
reproducing the endpoint behavior /? ~ 0. A reasonable cure for this problem (which we 
have not implemented here) is to add another basis function such as (1 — x^)^. 

Combining with T^(x) = nlln-i, and the orthogonality of the U^s under J dx{l — 
x^)^/^, we find that in this basis the 't Hooft integral has diagonal matrix elements 7r^n/2. 
The inner product is 



dx{l - x^)'U^-Y{x)lJn-\{x) 

K 

d9sin9sm.m9sva.n9 (A. 3) 

1 1 



l — {n—my^ 1 — (n+7Ti) 



n-m odd. 

The resulting generalized eigenvalue problem is very well behaved numerically and a 
basis size of 100 gives truncation errors less than 10~^° for the low lying masses. 
Changing the quark mass requires the matrix element 



< m\- -\n > = dxUm-i{x)Un-i{x) 

1 — x^ J 



sm mt' sm Tit' f \ a\ 

sin.9 l^-^J 



m — l 



2E 



n + m — 1 — 2A; 

We can adapt this to the higher kernels simply by computing their matrix elements. 
For an E^ theory we would need 

i<t>^., H,M ='/ (1^ tag (H^fjl^) (1 - -Y"U„.-,{.)(1 - yY"U„-,(.y)- 

(A.5) 
The integrand now has only a single pole and the singular cutoff goes smoothly to principal 
part evaluation as A ^ 0. This integral is particularly easy as we can write the log as the 
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sum of two terms, and in (say) the log(l — x)/(l + x) term do the y integral first, using 
(A.2) and taking the derivative. This wiU give 



{(j)^, Hs(j)n) = iT^n / dx{l -x^)^/^t/^_i(x)log ( , . I Un-i{x) - {n ^ m) 

= 27rm / (i^ sin m^ sin n^ log cot (n <-» m) 



Using GR 4.384.7 we find 



/ (i6' cos m6' log cot - = < I"*! (A. 7) 

Jo 2 10m even 

which applied to (|A.6| ) (keeping careful track of the absolute value sign) gives zero for 
m — n even (as required by parity) and for m — n odd gives 



2 I n — m n ~ m\ r, 2- , min(?7,, m) 



{(f)m, H3(j)n) = iT^ ] f ' = 27r^z sgu {n - m) -^^- (A.8) 

\\n — m\ n + mj n + m 

It is worth noting that the i in the bound state equation can be removed by a unitary 
transformation (which symmetrizes (|A.8| )). Formally this could be done in ( |3.24| ) but it 
somewhat obscures the structure. 
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